HLA-B evolutionary divergence is associated with outcomes after SARS-CoV-2 infection.

We examined the correlation between class I HLA evolutionary divergence (HED), a surrogate for the capacity to present different peptides, and the outcomes of 234 adult inpatients with confirmed SARS-CoV-2 infection. Genomic DNA was extracted from peripheral blood and genotyped by next-generation sequencing (NGS). HED scores for HLA class I (HLA-A, -B, and -C) genotypes were calculated using Grantham's distance. Higher HED scores for HLA-B, but not HLA-A or -C, are significantly associated with a decreased probability of poor outcomes including ICU admission, mechanical ventilation, and death (OR = 0.93; P = 0.04) in the univariate analysis. In the multivariate analysis, increased HLA-B HED score, younger age, and no comorbidity were independently associated with favorable outcomes (P = 0.02, P = 0.01, and P = 0.05, respectively). This finding is consistent with the notion that broader peptide repertoires presented by class I HLA may be beneficial in infection control.

HLA-Bw4 in association with KIR3DL1 favors natural killer cell-mediated protection against severe COVID-19.

Replicating SARS-CoV-2 has been shown to degrade HLA class I on target cells to evade the cytotoxic T-cell (CTL) response. HLA-I downregulation can be sensed by NK cells to unleash killer cell immunoglobulin-like receptor (KIR)-mediated self-inhibition by the cognate HLA-I ligands. Here, we investigated the impact of HLA and KIR genotypes and HLA-KIR combinations on COVID-19 outcome. We found that the peptide affinities of HLA alleles were not correlated with COVID-19 severity. The predicted poor binders for SARS-CoV-2 peptides belong to HLA-B subtypes that encode KIR ligands, including Bw4 and C1 (introduced by B*46:01), which have a small F pocket and cannot accommodate SARS-CoV-2 CTL epitopes. However, HLA-Bw4 weak binders were beneficial for COVID-19 outcome, and individuals lacking the HLA-Bw4 motif were at higher risk for serious illness from COVID-19. The presence of the HLA-Bw4 and KIR3DL1 combination had a 58.8% lower risk of developing severe COVID-19 (OR = 0.412, 95% CI = 0.187-0.904, p = 0.02). This suggests that HLA-Bw4 alleles that impair their ability to load SARS-CoV-2 peptides will become targets for NK-mediated destruction. Thus, we proposed that the synergistic responsiveness of CTLs and NK cells can efficiently control SARS-CoV-2 infection and replication, and NK-cell-mediated anti-SARS-CoV-2 immune responses being mostly involved in severe infection when the level of ORF8 is high enough to degrade HLA-I. The HLA-Bw4/KIR3DL1 genotype may be particularly important for East Asians undergoing COVID-19 who are enriched in HLA-Bw4-inhibitory KIR interactions and carry a high frequency of HLA-Bw4 alleles that bind poorly to coronavirus peptides.

HLA variants have different preferences to present proteins with specific molecular functions which are complemented in frequent haplotypes.

Human leukocyte antigen (HLA) genes are the most polymorphic loci in the human genome and code for proteins that play a key role in guiding adaptive immune responses by presenting foreign and self peptides (ligands) to T cells. Each person carries up to 6 HLA class I variants (maternal and paternal copies of HLA-A, HLA-B and HLA-C genes) and also multiple HLA class II variants, which cumulatively define the landscape of peptides presented to T cells. Each HLA variant has its own repertoire of presented peptides with a certain sequence motif which is mainly defined by peptide anchor residues (typically the second and the last positions for HLA class I ligands) forming key interactions with the peptide-binding groove of HLA. In this study, we aimed to characterize HLA binding preferences in terms of molecular functions of presented proteins. To focus on the ligand presentation bias introduced specifically by HLA-peptide interaction we performed large-scale in silico predictions of binding of all peptides from human proteome for a wide range of HLA variants and established which functions are characteristic for proteins that are more or less preferentially presented by different HLA variants using statistical calculations and gene ontology (GO) analysis. We demonstrated marked distinctions between HLA variants in molecular functions of preferentially presented proteins (e.g. some HLA variants preferentially present membrane and receptor proteins, while others - ribosomal and DNA-binding proteins) and reduced presentation of extracellular matrix and collagen proteins by the majority of HLA variants. To explain these observations we demonstrated that HLA preferentially presents proteins enriched in amino acids which are required as anchor residues for the particular HLA variant. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Taking into consideration that HLA alleles are inherited in haplotypes, we hypothesized that haplotypes composed of a combination of HLA variants with different presentation preferences should be more advantageous as they allow presenting a larger repertoire of peptides and avoiding holes in immunopeptidome. Indeed, we demonstrated that HLA-A/HLA-B and HLA-A/HLA-C haplotypes which have a high frequency in the human population are comprised of HLA variants that are more distinct in terms of functions of preferentially presented proteins than the control pairs.

Distribution of KIR Genes and Their HLA Ligands in Different Viral Infectious Diseases: Frequency Study in Sicilian Population.

Natural killer (NK) cells play a role in defence against viral infections by killing infected cells or by producing cytokines and interacting with adaptive immune cells. Killer immunoglobulin-like receptors (KIRs) regulate the activation of NK cells through their interaction with human leucocyte antigens (HLA). Ninety-six Sicilian patients positive to Human Immunodeficiency Virus-1 (HIV) and ninety-two Sicilian patients positive to SARS-CoV-2 were genotyped for KIRs and their HLA ligands. We also included fifty-six Sicilian patients with chronic hepatitis B (CHB) already recruited in our previous study. The aim of this study was to compare the distribution of KIR-HLA genes/groups of these three different infected populations with healthy Sicilian donors from the literature. We showed that the inhibitory KIR3DL1 gene and the KIR3DL1/HLA-B Bw4 pairing were more prevalent in individual CHB. At the same time, the frequency of HLA-C2 was increased in CHB compared to other groups. In contrast, the HLA-C1 ligand seems to have no contribution to CHB progression whereas it was significantly higher in COVID-19 and HIV-positive than healthy controls. These results suggest that specific KIR-HLA combinations can predict the outcome/susceptibility of these viral infections and allows to plan successful customized therapeutic strategies.

HLA Allele Frequencies and Association with Severity of COVID-19 Infection in Northern Italian Patients.

HLA allelic distribution was analysed in a cohort of 96 Northern Italian subjects (53M/43F) (mean age 59.9 ± 13.3 years) from Lombardy who developed COVID-19 during the first two pandemic waves to investigate possible correlations between HLA molecules and disease severity. An important role of HLA- B and HLA-C loci in modulating the clinical severity of COVID-19 disease was identified. In particular, the HLA-B07 supertype was observed to be associated with a significant risk for severe disease; conversely, the HLA-B27 supertype and C*12:02 allele played a protective role as they were associated with milder disease. These associations were confirmed after applying a multinomial regression analysis to adjust the correlation for age, gender and comorbidities with COVID-19 severity. Though the power of results is limited by the small sample size, data herein contribute to shedding light on the role played by genetic background in COVID-19 infection.

Discovery of a novel HLA-B*40 allele, HLA-B*40:02:01:30 in a Japanese individual.

HLA-B*40:02:01:30 differs from HLA-B*40:02:01:01 by one nucleotide in 3'UTR at position 2869.

Distribution of HLA-A, -B, -C, -DRB1, -DQB1, -DPB1 allele frequencies in patients with COVID-19 bilateral pneumonia in Russians, living in the Chelyabinsk region (Russia).

In this population-based case-control study conducted in the Chelyabinsk region of Russia, we examined the distribution of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1, in a group of 100 patients with confirmed COVID-19 bilateral pneumonia. Typing was performed by NGS and statistical calculations were carried out with the Arlequin program. HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 alleles were compared between patients with COVID-19 and 99 healthy controls. We identified that COVID-19 susceptibility is associated with alleles and genotypes rs9277534A (disequilibrium with HLA-DPB1*02:01, -02:02, -04:01, -04:02, -17:01 alleles) with low expression of protein products HLA-DPB1 (pc < 0.028) and homozygosity at HLA-C*04 (p = 0.024, pc = 0.312). Allele HLA-A*01:01 was decreased in a group of patients with severe forms of bilateral pneumonia, and therefore it may be considered as a protective factor for the development of severe symptoms of COVID-19 (p = 0.009, pc = 0.225). Our studies provide further evidence for the functional association between HLA genes and COVID-19.

Association of Human Leukocyte Antigen Genotypes with Severe Acute Respiratory Syndrome Coronavirus 2 Vaccine-Induced Subacute Thyroiditis.

Background: Despite mass vaccination, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine-induced subacute thyroiditis (SAT) is rarely seen as a complication. The reason why some individuals are susceptible to developing vaccine-induced SAT is not known. SAT develops in genetically predisposed individuals who carry specific human leukocyte antigen (HLA) haplotypes. It is unknown whether specific HLA alleles are associated with SARS-CoV-2 vaccine-induced SAT. Objective: This study compared the HLA profiles of patients with SARS-CoV-2 vaccine-induced SAT to controls, to assess whether there is an association between specific HLA genotypes and development of SAT. The relationship between HLA genotypes and the clinical course of SARS-CoV-2 vaccine-induced SAT was also evaluated. Methods: A case-control study was conducted in a Turkish tertiary care center. Fourteen patients with SARS-CoV-2 vaccine-induced SAT and 100 healthy controls were included. HLA-A, HLA-B, HLA-C, HLA-DQB1, and HLA-DRB1 frequencies were analyzed by next-generation sequencing. Results: The frequencies of HLA-B*35 and HLA-C*04 alleles were significantly higher in SARS-CoV-2 vaccine-induced SAT cohort when compared with controls (HLA-B*35: 13 [93%] vs. 40 [40%], p < 0.001; HLA-C*04: 13 [93%] vs. 43 [43%], p < 0.001, respectively). More severe thyrotoxicosis was seen in patients having HLA-B*35 and HLA-C*04 homozygous alleles (free thyroxine: 4.47 ng/dL [3.77-5.18] vs. 1.41 ng/dL [1.22-2.63], p = 0.048). Inflammation tended to be more severe in homozygous patients (C-reactive protein: 28.2 mg/dL [13.6-42.9] vs. 4.8 [1.2-10.5], p = 0.07). Conclusions: The frequencies of HLA-B*35 and HLA-C*04 alleles were higher in SARS-CoV-2 vaccine-induced SAT compared with controls. Homozygosity for HLA-B*35 and HLA-C*04 was associated with thyrotoxicosis and a greater inflammatory reaction. Our findings should be confirmed in studies of other populations.

Allelic imbalance of HLA-B expression in human lung cells infected with coronavirus and other respiratory viruses.

The human leucocyte antigen (HLA) loci have been widely characterized to be associated with viral infectious diseases using either HLA allele frequency-based association or in silico predicted studies. However, there is less experimental evidence to link the HLA alleles with COVID-19 and other respiratory infectious diseases, particularly in the lung cells. To examine the role of HLA alleles in response to coronavirus and other respiratory viral infections in disease-relevant cells, we designed a two-stage study by integrating publicly accessible RNA-seq data sets, and performed allelic expression (AE) analysis on heterozygous HLA genotypes. We discovered an increased AE pattern accompanied with overexpression of HLA-B gene in SARS-CoV-2-infected human lung epithelial cells. Analysis of independent data sets verified the respiratory virus-induced AE of HLA-B gene in lung cells and tissues. The results were further experimentally validated in cultured lung cells infected with SARS-CoV-2. We further uncovered that the antiviral cytokine IFNß contribute to AE of the HLA-B gene in lung cells. Our analyses provide a new insight into allelic influence on the HLA expression in association with SARS-CoV-2 and other common viral infectious diseases.

The Association between Proportion HLA-BW4 or HLA-BW6 May Causes Immunity Failure in COVID-19.

No Abstract.

Association between the HLA genotype and the severity of COVID-19 infection among South Asians.

Regional variations are found in the incidence and severity of the COVID-19 infection. Human leukocyte antigen (HLA) polymorphism is one of the genetic factors that might have an impact on the outcome of the disease. This study explored the association between the HLA genotype and the severity of COVID-19 among patients from South Asia. Blood samples from 95 Asians (Bangladeshis, Indians, and Pakistanis) with COVID-19 were collected. The patients were divided according to the severity of their infection: mild (N = 64), severe (N = 31), and fatal (N = 20). DNA was extracted from all samples, and HLA genotyping was performed for both class I (A, B, and C) and class II (DRB1, DQA1, and DQB1) using the PCR-rSSO (polymerase chain reaction-reverse sequence-specific oligonucleotide) molecular method. The frequency of HLA-B*51 was significantly higher among patients in the fatal group than among those in the mild infection group (15% vs. 4.7%, p = 0.027). Additionally, the frequency of HLA-B*35 was significantly higher in the mild infection group than in the fatal group (21.1% vs. 7.5%, p = 0.050 [a borderline p-value]). In terms of HLA class II, DRB1*13 was significantly observed in the fatal group than in the mild infection group (17.5% vs. 11.3%, p = 0.049). However, the p-value for all alleles became insignificant after a statistical correction for the p-values (pc = 0.216, pc = 0.4, and pc = 0.49, respectively). Compared with all published data, this study highlights that the association between the HLA system and the COVID-19 outcome might be ethnic-dependent. Genetic variation between populations must be examined on a wider scale to assess infection prognosis and vaccine effectiveness.

Predictive immunogenetic markers in COVID-19.

Immunorelevant genes are among the most probable modulators of coronavirus disease 2019 (COVID-19) progression and prognosis. However, in the few months of the pandemic, data generated on host genetics has been scarce. The present study retrieved data sets of HLA-B alleles, KIR genes and functional single nucleotide polymorphisms (SNPs) in cytokines related to COVID-19 cytokine storm from two publicly available databases: Allele Frequency Net Database and Ensembl, and correlated these frequency data with Case Fatality Rate (CFR) and Daily Death Rates (DDR) across countries. Correlations of eight HLA-B alleles and polymorphisms in three cytokine genes (IL6, IL10, and IL12B) were observed and were mainly associated with DDR. Additionally, HLA-B correlations suggest that differences in allele affinities to SARS-CoV-2 peptides are also associated with DDR. These results may provide rationale for future host genetic marker surveys on COVID-19.

Class I HLA Allele Predicted Restricted Antigenic Coverages for Spike and Nucleocapsid Proteins Are Associated With Deaths Related to COVID-19.

The world is dealing with one of the worst pandemics ever. SARS-CoV-2 is the etiological agent of COVID-19 that has already spread to more than 200 countries. However, infectivity, severity, and mortality rates do not affect all countries equally. Here we consider 140 HLA alleles and extensively investigate the landscape of 3,723 potential HLA-I A and B restricted SARS-CoV-2-derived antigens and how 37 countries in the world are predicted to respond to those peptides considering their HLA-I distribution frequencies. The clustering of HLA-A and HLA-B allele frequencies partially separates most countries with the lowest number of deaths per million inhabitants from the other countries. We further correlated the patterns of in silico predicted population coverage and epidemiological data. The number of deaths per million inhabitants correlates to the predicted antigen coverage of S and N derived peptides and its module is influenced if a given set of frequent or rare HLA alleles are analyzed in a given population. Moreover, we highlighted a potential risk group carrying HLAs associated with an elevated number of deaths per million inhabitants. In addition, we identified three potential antigens bearing at least one amino acid of the four-length insertion that differentiates SARS-CoV-2 from previous coronavirus strains. We believe these data can contribute to the search for peptides with the potential to be used in vaccine strategies considering the role of herd immunity to hamper the spread of the disease. Importantly, to the best of our knowledge, this work is the first to use a populational approach in association with COVID-19 outcome.

Association of HLA-B22 serotype with SARS-CoV-2 susceptibility in Hong Kong Chinese patients.

The coronavirus disease 2019 (COVID-19) is a highly infectious disease caused by SARS-CoV-2. Since its first report in December 2019, COVID-19 has evolved into a global pandemic causing massive healthcare and socioeconomic challenges. HLA system is critical in mediating anti-viral immunity and recent studies have suggested preferential involvement of HLA-B in COVID-19 susceptibility. Here, by investigating the HLA-B genotypes in 190 unrelated Chinese patients with confirmed COVID-19, we identified a significant positive association between the B22 serotype and SARS-CoV-2 infection (p = 0.002, Bonferroni-corrected p = 0.032). Notably, the B22 serotype has been consistently linked to susceptibility to other viral infections. These data not only shed new insights into SARS-CoV-2 pathogenesis and vaccine development but also guide better infection prevention/control.

Correlation of the two most frequent HLA haplotypes in the Italian population to the differential regional incidence of Covid-19.

BACKGROUND: Understanding how HLA polymorphisms may affect both susceptibility, course and severity of Covid-19 infection could help both at the clinical level to identify individuals at higher risk from the disease and at the epidemiological one to explain the differences in the epidemic trend among countries or even within a specific country. Covid-19 disease in Italy showed a peculiar geographical distribution from the northern most affected regions to the southern ones only slightly touched. METHODS: In this study we analysed the regional frequencies for the most common Italian haplotypes from the Italian Bone Marrow Donor Registry (HLA-A, -B, -C and -DRB1 at four-digit level). Then we performed Pearson correlation analyses among regional haplotypes estimated frequency in the population and Covid-19 incidence and mortality. RESULTS: In this study we found that the two most frequent HLA haplotypes in the Italian population, HLA-A*:01:01g-B*08:01 g-C*07:01g-DRB1*03:01g and HLA-A*02.01g-B*18.01g-C*07.01g-DRB1*11.04g, had a regional distribution overlapping that of Covid-19 and showed respectively a positive (suggestive of susceptibility) and negative (suggestive of protection) significant correlation with both Covid-19 incidence and mortality. CONCLUSIONS: Based on these results, in order to define such HLA haplotypes as a factor effectively associated to the disease susceptibility, the creation of national networks that can collect patients' samples from all regions for HLA typing should be highly encouraged.

A role for human leucocyte antigens in the susceptibility to SARS-Cov-2 infection observed in transplant patients.

We analysed data from 80 patients who tested positive for SARS-CoV-2 RNA who had previously been HLA typed to support transplantation. Data were combined from two adjacent centres in Manchester and Leeds to achieve a sufficient number for early analysis. HLA frequencies observed were compared against two control populations: first, against published frequencies in a UK deceased donor population (n = 10,000) representing the target population of the virus, and second, using a cohort of individuals from the combined transplant waiting lists of both centres (n = 308), representing a comparator group of unaffected individuals of the same demographic. We report a significant HLA association with HLA- DQB1*06 (53% vs. 36%; p < .012; OR 1.96; 95% CI 1.94-3.22) and infection. A bias towards an increased representation of HLA-A*26, HLA-DRB1*15, HLA-DRB1*10 and DRB1*11 was also noted but these were either only significant using the UK donor controls, or did not remain significant after correction for multiple tests. Likewise, HLA-A*02, HLA-B*44 and HLA-C*05 may exert a protective effect, but these associations did not remain significant after correction for multiple tests. This is relevant information for the clinical management of patients in the setting of the current SARS-CoV-2 pandemic and potentially in risk-assessing staff interactions with infected patients.

SARS-CoV-2 genomic variations associated with mortality rate of COVID-19.

The coronavirus disease 2019 (COVID-19) outbreak, caused by SARS-CoV-2, has rapidly expanded to a global pandemic. However, numbers of infected cases, deaths, and mortality rates related to COVID-19 vary from country to country. Although many studies were conducted, the reasons of these differences have not been clarified. In this study, we comprehensively investigated 12,343 SARS-CoV-2 genome sequences isolated from patients/individuals in six geographic areas and identified a total of 1234 mutations by comparing with the reference SARS-CoV-2 sequence. Through a hierarchical clustering based on the mutant frequencies, we classified the 28 countries into three clusters showing different fatality rates of COVID-19. In correlation analyses, we identified that ORF1ab 4715L and S protein 614G variants, which are in a strong linkage disequilibrium, showed significant positive correlations with fatality rates (r = 0.41, P = 0.029 and r = 0.43, P = 0.022, respectively). We found that BCG-vaccination status significantly associated with the fatality rates as well as number of infected cases. In BCG-vaccinated countries, the frequency of the S 614G variant had a trend of association with the higher fatality rate. We also found that the frequency of several HLA alleles, including HLA-A*11:01, were significantly associated with the fatality rates, although these factors were associated with number of infected cases and not an independent factor to affect fatality rate in each country. Our findings suggest that SARS-CoV-2 mutations as well as BCG-vaccination status and a host genetic factor, HLA genotypes might affect the susceptibility to SARS-CoV-2 infection or severity of COVID-19.